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Abstract 

We update the electroweak study of the predictions of the Minimal Supersymmetric 
Standard Model (MSSM) including the recent results on the muon anomalous magnetic 
moment, the weak boson masses, and the final precision data on the Z boson param- 
eters from LEP and SLC. We find that the region of the parameter space where the 
slepton masses are a few hundred GeV is favored from the muon g — 2 for tan (3 ^ 10, 
whereas for tan/3 ~ 50 heavier slepton mass up to ~ 1000 GeV can account for the 
reported 3.2 a difference between its experimental value and the Standard Model (SM) 
prediction. As for the electroweak measurements, the SM gives a good description, and 
the sfermions lighter than 200 GeV tend to make the fit worse. We find, however, that 
sleptons as light as 100 to 200 GeV are favored also from the electroweak data, if we 
leave out the jet asymmetry data that do not agree with the leptonic asymmetry data. 
We extend the survey of the preferred MSSM parameters by including the constraints 
from the 6—7-57 transition, and find favorable scenarios in the minimal supergravity, 
gauge-, and mirage-mediation models of supersymmetry breaking. 



1 Introduction 



Despite anticipation that physics beyond the Standard Model (SM) should show up at ener- 
gies just above the reach of the present collider experiments, we have so far been unsuccessful 
in identifying the nature of new physics. Supersymmetric (SUSY) extensions of the SM, in 
particular, the minimal supersymmetric standard model (MSSM) has several attractive fea- 
tures, such as the unification of the three gauge couplings. In Ref. P, two of us performed 
a comprehensive study of constraints on the MSSM parameters from precision electroweak 
(EW) experiments by using the data published in 1999 [2]. It was found that almost all 
SUSY particle masses are constrained to be larger than a few 100 GeV, except for the light 
wino-like chargino whose contribution improved the SM fit slightly. 

Since then, there have been several improvements in the EW measurements and theo- 
retical analyses. Most notably, the LEP results have been finalized [3l H] and the estimate 
of the running QED coupling constant at the Z boson mass scale has been improved by the 
contribution from the BES experiment O [6] . Also, the measurements of the muon anoma- 
lous magnetic moment at BNL have been finalized [7], and a 3.2 a discrepancy from the 
SM prediction has been reported |H] (see Refs. [9-11] for the determination of the hadronic 
contribution to the muon g — 2 from other groups). Those changes in the EW data, the W 
boson mass data [121 [13] ^iiid the top-quark mass data [131 E] l^ad to new constraints on 
the MSSM parameters. On theory side, a number of groups have studied the EW fits in the 
SM and in the MSSM: for reviews, see e. g. Ref. [U]. Refs. [16-18] update the constraints 
on the Higgs boson mass in the SM using the EW data as well as the LEP and the Teva- 
tron data for direct searches for the SM Higgs boson. Ref. [19] provides a comprehensive 
analysis of the MSSM EW precision fits using the state-of-the-art multi-loop calculations 
of the EW precision observables [20] . The EW precision fits in the minimal supergravity 
mediated SUSY breaking (mSUGRA) scenarios have been studied in many papers [21-47]. 
They are also studied in the split SUSY model [121 SH], in the gauge- [271 [il ESI US SZl SS] , 
anomaly- [2H1 ESI 112] , moduli- [2Sl[Sn], and radion mediations [12], in the supergravity mod- 
els with non- universal gaugino masses [51] and with non-universal Higgs masses [261 152] and 
in the 25-parameter "phenomenological MSSM" model [53]. The EW observables in the 
MSSM are also studied in attempts to accommodate the NuTeV anomaly [SH |55] and the 
discrepancy between the effective Weinberg angles extracted from the leptonic and hadronic 
asymmetries [SS]. 

Another advance on the theory side is a deeper understanding of the jet asymmetries 
measured at the LEP experiments. In Ref. [57], QCD radiative corrections to the jet asym- 
metries are studied, and it is argued that in the final report of the LEP experiments [3] the 
systematic errors of the jet asymmetries might have been underestimated. This is also a 
part of the motivation to revisit the EW fits in this paper. 

In this paper, we present quantitative results based on the muon g — 1 and the final EW 
data from LEP and SLC [3lll]. In section 2, we discuss the MSSM contribution to the muon 
(7 — 2, and identify its preferred parameter region. In section 3, we give our parametrization 
of the EW observables at the Z-pole and the mass and the width of the W boson in the SM. 
In section 4, we briefly review the MSSM contributions to the EW precision observables. In 
section 5, we explore a few SUSY breaking models and identify several preferred scenarios 
that can accommodate also the 6 — ?■ S7 rate. Section 6 gives summary and discussions. 
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2 The muon ^ - 2 vs MSSM 



2.1 The muon g - 2 in the MSSM 



The muon g — 2 is a precisely measured quantity, and hence it is an excellent probe of new 
physics at the TeV scale. The measurement at BNL was finalized in 2006 [7j. After including 
a small shift in the value of the proton-to-muon magnetic ratio reported since then [58] , the 
experimental value is 



= 11659208.9(6.3) x 10 



-10 



As for the SM prediction, the recent improvement in the low-energy e~^e — )■ hadrons data 
from BaBar [SH], BES H, CMD-2 [BU], KLOE [SQ E2] and SND P] allows us to reduce the 
uncertainty. In Ref. [8], the SM prediction has been evaluated as (after correcting a typo in 
Ref. [8]), 



11659183.0(5.1) X 10^ 



(2) 



which includes all the major updates on the e~^e~ — )■ hadrons data, and adopts the esti- 
mate |6l], a^^^^'^ = 10.5(2.6) x 10~^°, for the light-by-light contribution. Other independent 
analyses [9l[Tni[lI] based on the e~^e~ data give similar estimates. Hence, the observed value 
of the muon g — 2 is larger than the SM prediction by 



Sa,. 



cxp 



(25.9 ±8.1) X 10^ 



-10 



(3) 



which differs from zero by 3.2 cr. It is tempting to interpret the difference as a contribution of 
new particles with the muon quantum number, such as smuons in the SUSY SM. Throughout 
this article we assume that the MSSM contribution accounts for this discrepancy. 

In the MSSM, the contribution to the muon g — 2 has been calculated up to and including 
the two-loop level [65]. In view of the smallness of the two- loop contribution, in this paper 
we restrict our analyses in the one-loop approximation^. 

At one-loop, the MSSM contribution comes from the chargino contribution a^(x~) and 
the neutralino contribution a^(x°). The relevant one-loop expressions in the notation of 
Ref. [1] are found e.g. in Refs. [6 




(4a) 



(4b) 



^ The effect from tlie tan /3-enlianced resummation in tfie muon-muon-Higgs vertex can cliange tlie MSSM 
contribution by about 10% [551, which does not affect the conclusions in the present paper significantly. 
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(b) 





Figure 1: The SUSY contributions to the muon g — 2 which give the leading terms of the 
expansion in m^/msusY- The photon (wavy hne) is attached to all the charged particles. 



where 



G2(X) 

GJx) 



Gi(x) = — ; — \(x — — 5x — 2) + Gxlnxl , 

IV ; i2(x- 1)4 Lv A ; j , 

: — ; \(x — l)(2x^ + 5x — 1) — Gx^lnxl 

12 X - 1 4 LV A ; J 



2{x - 1)3 
1 

^2rx- 1)3 



[(x- l)(x-3) + 21nx] , 
[{x — l)(x + 1) — 2xlna;] . 



(5a) 
(5b) 
(5c) 
(5d) 



Even though these expressions are useful for numerical calculations, they are not particularly 
illuminating for the purpose of understanding their dependences on the SUSY parameters. 
The main disadvantage of the above expressions is that they are written in terms of the 
mass eigenstates, in terms of which the dependences on the SUSY breaking parameters are 
hidden by the electroweak symmetry breaking that causes complex mixings. 

In the weak eigenstates, the structure of the one-loop contributions becomes much more 
transparent. This simplification occurs since the expressions in the weak eigenstates are 
equivalent to the m^/msusY expansion, where msusY is the typical SUSY breaking mass 
scale. The price we have to pay is that the leading terms in the expansion are not useful 
when msusY ~ f^z- However, we will find below that this expansion is very useful when 
analyzing the SUSY parameter dependence. 

The leading terms in the mz /msijsY expansion are given by the five diagrams (a) to (e) 
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Figure 2: The behaviors of the functions Fa{x,y) and Ff,{x,y), which appear in SUSY 
contributions to the muon g — 2, for y = x,2x,3x. 



in Fig. [H whose contributions can be expressed compactly as 

a^{W-H, v^) = — -J Fa (^^, ^ J , (6a) 

a,{Br^,-~^^)= ^^^^^^i^J^,^ , (6b) 



a,{B-H,~^L)= FJ^,^], (6c) 
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a,(M,/i«) = -15^^45 Fj-^,:^], (6e 
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respectively. The functions Fa{x,y) and Fh{x,y) are defined as: 

P / N _ (^3(3:) -G3(y) p , ^ _ <:?4(a:) - ^4(1/) 

^a(a;,?/) = , Ffe(x,?/) = , (7) 

X — y X — y 

which are symmetric under exchange of the two arguments. The functions Gs{x) and 6*4 (x) 
are monotonically decreasing for < x < 00, and hence the functions Fa{x, y) and Fh{x, y) 
are positive for all positive x and y. In Fig. [21 we show the behaviors of Fa^b{x,nx) = 
Fafi{nx, x) for 77, = 1, . . . , 3. 

The expressions (l6^)-([6^) allow us to make a few general observations on the SUSY 
parameter dependences. The first one is that the contributions from the diagrams (a)-(c) 
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in Fig. [U are positive, while those from the diagrams (d) and (e) are negative for M2/i and 
Mi/_i > 0. In addition, if the mass-sphtting in the (i>^,/ii) doublet is small we can conclude 
that the sum of the diagrams (a) and (d), or that of Eqs. ([6^) and ([6li) is always positive, 
because Fa{x,y) is always larger than Fb{x,y) for the same arguments as shown in Fig. [2l 
If the contributions from the diagrams (b), (c), and (e) are suppressed, being proportional 
to gy = g^taiL^Ow, we can conclude that M2fi should be positive in order to explain the 
positive deviation of the muon g — 2 from the SM prediction in Eq. ([3]). 

In fact, the only way to obtain a positive MSSM contribution for Mi/i < and M2/i < 
is to make the contribution from the diagram (e), Eq. ([6^), dominates over the others, which 

requires (M2/m|)Fa(M|/m|, /iV^|) < tan^ 6'vi/(Mi/m|)F6(Mf/m|, /iV^|)- This situa- 
tion is realized for \Mi\,m^ -C \M2\,mi: for instance, when (tan/3, /i, m^,, mj^, M2, Mi) = 
(50, —300, 120, 1200, 1000, 100), where the mass dimension is counted in units of GeV, the 
diagram (e) dominates over the others, and the predicted value of is within the 1-a favored 
region. In this example, an Mi/i < solution is realized by setting M2/M1 = m^/m^ = 10. 
We find that all the solutions with Mi/i < require both IM2/M1I ^ 1 and rrii/m^ ^ 1. 
However, as discussed later, scenarios with heavy left-handed sleptons do not lead to a 
significant signal in the electroweak precision measurements, and we do not consider such 
scenarios hereafter. 

Another interesting observation is that in the case where /i is large as compared to the 
SUSY breaking soft mass parameters, the contributions from the diagrams (a) and (c)-(e) 
are suppressed by l//i, as manifestly shown by the Higgsino propagator in the diagrams 
Figs. [T] (a) and (c)-(e). On the other hand, the diagram (b) is proportional to fi, which 
makes this contribution important. 

The above observations can be made explicit by neglecting the D-term and the F-term 
contributions to the slepton mass matrices, which makes m^^ = nip = mi and m^^ = m^, 
where mi and m^ are the left-handed and the right-handed slepton soft SUSY breaking 
masses, respectively. Then, Eqs. (|^ -(lBel) simplify as 




(8a) 




(8b) 



(8c) 




(8d) 



ref 



tan^ 9' 




) 



(8e) 



where 



a 



ref _ 

'P- = 



g^ m'j^M2fita.n (3 



(9) 
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Only for the above specific scenario witfi mi/m^ ^ 1 and IM2/M1I 3> 1, Eq. ( IHel) dominates 
over the others and Mi/i < is required to account for the positive shift, Eq. ([3]). In the very 
heavy Higgsino case with |/i/Mi|, |/i/M2| ^ 1, Eq. ( l8b|) dominates, and Mi/i > is necessary 
to account for the positive shift. Except for the above two cases, Eq. ( !8a|l dominates, and 
M2IJ, > follows. 

We note that the above discussion on the favored signs of Mifi and M2fi are valid equally 
for both positive and negative Mi and M2 since the signs of Mi and M2 enter only through 
the combinations Min and M2/i. In particular, we can use the same argument to identify 
the favored signs of Mi/i and M2/i also for a model with sgn(Mi) = — sgn(M2), which is 
realized in some parameter regions of the mixed moduli and anomaly mediation model |69] . 

A comment on the SUSY breaking tri-linear coupling is in order here. In the one- 
loop order, the parameter enters only through the jli-fj'R mixing in the combination 
(A* — /itan/3) in the diagram Fig. [T] (b). In order for to affect the muon g — 2 sig- 
nificantly. Fig. [1] (b) should dominate over the others, but it implies large \n\ according to 
the above discussions. Since there is no attractive SUSY breaking scenarios that lead to a 
huge magnitude of 1^4^ | ~ tan/3, and since moderate values of |y4^| ^ do not affect the 
qualitative behavior, we set A^ = in the following, unless otherwise stateci^. 

In Fig. [3] we show contour plots of the SUSY contribution to the muon — 2 as a function 
of the right-handed slepton mass and the SU{2)l gaugino mass M2. In the figures, 
we examine the four combinations of tan/3 and fi; tan/3 = 10 or 50 and fi = 200 or 800 
GeV. As for the other relevant SUSY parameters, we assume the gaugino mass unification 
condition, M1/M2 = (5/3) tan^ 6w, and the left-handed slepton mass is taken to be equal to 
the right-handed slepton mass, for simplicity. Under these assumptions, we plot, from the 
lower left corner, the +3a, +2a, +la, —la and —2a contours. (We do not plot —3a curves 
since those lie outside of the panel for Fig. |3]^a), (b), (c) and (d)). The area surrounded by 
the ±1(T contours is shaded. The vertical dotted lines give the lower bound on M2 beyond 
which the lightest chargino mass becomes less than the experimental limit, 103.5 GeV [7Tj . 
Also shown by the horizontal dotted lines are the stau mass limit, 81.9 GeV [13]. We also 
show the contour of a constant "electroweak factor" , defined as the squared sum of the 
SUSY contributions to the electroweak Sz- and T^-parameters [1] and to mw] 

. 2 f {^Sz)svsY y , /^ (AT^)susy V , f {Amw)s\JSY Y 
^^^^=[ 0.04 ) +1, 0.04 ) +1 30MeV ) ' ^'^^ 

The parameter regions in the left-bottom side of the dash-dotted curves give ^Xew > 0.5, 
and there can be a hint of SUSY contributions to the electroweak observables in the near 
futur^. See discussions in Section 4 for details. 

In the figures, the sample points MSSM1-MSSM6 denoted by the cross symbols (x), are 
chosen in the 1-a allowed region of the muon g — 2 where Axew > 0.5. They represent the 
two extreme cases where either a slepton (2,4,6) or an ino (1,3,5) is light. Our choices of fi and 
M2 allow us to cover the cases where the lighter chargino is wino-like (3,4,5), Higgsino-like 

^ For the possibility of explaining the muon g — 2 anomaly using a large A term, see e.g. Ref. [70] . 

3 The "dips" in Figs. 3 (a) and (c) at M2 ~ 120GeV and 350 GeV ^ ^ 400_GcV happen since around 
this region there is a destructive interference between the combination "AT — ASq / a" and the ARz term 
in ATz (see eq. (|16bp ). which makes the contribution from ATz to Ax|^ small. 
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Figure 3: The muon g — 2, plotted against M2 (the SU{2)l gaugino mass) and (the right- 
handed smuon soft SUSY breaking mass) for tan /3 = 10 (top two panels) and tan /3 = 50 
(bottom two panels), and for /i = 200 GeV (left two panels) and fi = 800 GeV (right two 
panels). The curves are, from the lower left corner, +3a, +2cr, +lcr, —la and —2a contour 
for the difference (5a^ between the data and the SM prediction. The region on the left- 
hand side of the vertical dotted line is excluded by the chargino mass limit m-- > 103.5 
GeV [7T], and the region below the horizontal dotted line is ruled out by the stau mass limit 
rrif^ > 81.9 GeV [13]. The region below or in the left-hand side of the dash-dotted curve 
gives Axew > 0-^ contribution to the electroweak observables, see Eq. ( ITOj) . The sample 
points discussed in the main text are marked by the crosses (x). In the figures, we assume 
A^ = 0,mi = m^ and M1/M2 = (5/3) tan^ Ow 
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(2), or mixed (1,6). We use these six MSSM sample points in the study of the electroweak 
observables in later sections. The choice that there is either a light ino or a light slepton is 
interesting in view of the electroweak study since if there is a light new particle with nonzero 
electroweak quantum numbers, the contribution to the Sz- and T^-parameters [1] and to 
mw can be significant. 

For these sample points, we present separate contributions from individual diagrams and 
their sum in Tabled! In the table, we show the sum of the five terms in the column '(a)-(e)', 
and the total SUSY contribution obtained without using the mz/msvsY expansion in the 
column under 'total'. By comparing these two columns, we confirm for all the six cases that 
the leading terms in the mz/msnsY expansion give an excellent approximation. The last 
column shows the pull factor 

puii^(^°^)7;-(^°^>. („) 

Ada^ 

where the mean (5a^) and the error A5a^ are taken from Eq. ([3]). 

For tan/3 = 10 and fi = 200 GeV shown in Fig. ^a), the 1-a allowed region is very 
roughly given by 550GeV< 1.7m^+M2 < 800 GeV. Here, the dia gram (a) gives the dominant 
contribution, as can be read off from the rows 1 and 2 in Tabled], respectively, for the MSSM 
points 1 and 2. 

For a larger fi like /i = 800 GeV shown in Fig. [3](b), the allowed region becomes narrower, 
partly due to the heavier Higgsinos, and partly because of the constraint from the stau mass 
lower bound > 81.9 GeV [13]. In the parameter region of Fig. |3]^b), the diagram (b) 
becomes more important because all the other diagrams are suppressed by the heavy Higgsino 
mass, as discussed above. This can be explicitly verified from Table dl in the rows 3 and 4. 
In the MSSM point 3, the suppression of the diagram (a) by l//i is not strong enough and 
hence the diagram (a) is still as important as (b). In the MSSM point 4, the diagram (a) is 
less important since it is suppressed not only by l//i but also by I/M2. 

In the tan /3 = 50 case, the allowed parameter space becomes much wider as shown in 
Fig. El^c) for /i = 200 GeV and Fig. El^d) for /x = 800GeV. When /x ~ 200 GeV, the favored 
SUSY masses are so large that there is no region which satisfies both the 1-a favored range 
of the muon g — 2 and Ax^w > O-^- When fi ~ 800 GeV, in Fig. M^d), there appear 
two distinct regions of the parameters that satisfy both conditions. In the small M2 region 
around the MSSM point 5, M2 ^ 100 GeV is allowed by the lighter chargino mass constraint, 
and M2<150 GeV gives AxIyv > I^ small region around the MSSM point 6, 
although 300 GeV, the lighter stau is as light as ~100 GeV and it gives a sizable 

contribution to the AT-parameter, which makes Ax^^v non-negligible. As for the muon 
g — 2/m. the MSSM point 5, the large slepton mass suppresses the diagram (b) despite large 
/i, and the light wino makes the diagram (a) dominate over the other contributions. In the 
MSSM point 6, the large /i enhancement of the diagram (b) is more effective, while the 
diagram (a) is suppressed by the large wino mass M2. As a result, the contributions from 
the diagrams (a) and (b) are comparable. 

In summary, for the scenarios in which the SUSY contribution to the muon g — 2 can 
be tested by the electroweak precision study, the diagrams (a) and/or (b) give dominant 
contribution to the muon g — 2. 
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No. 


tan (3 


/i 


Ms 




(a) 


(b) 


(c) 


(d) 


(e) 


(a)-(e) 


total 


pull 


1 


10 


200 


150 


300 


29.6 


1.1 


0.7 


-2.9 


-1.3 


27.2 


25.0 


-0.1 


2 


10 


200 


450 


120 


27.5 


8.8 


3.3 


-7.1 


-6.7 


25.9 


25.9 


0.0 


3 


10 


800 


150 


200 


14.3 


16.2 


0.6 


-2.7 


-1.3 


27.1 


27.1 


0.1 


4 


10 


800 


500 


150 


6.9 


21.3 


1.0 


-2.5 


-2.1 


24.7 


24.3 


-0.2 


5 


50 


800 


150 


550 


26.9 


2.4 


0.5 


-2.6 


-1.0 


26.3 


26.0 


0.0 


6 


50 


800 


900 


280 


18.0 


18.0 


2.5 


-5.9 


-5.1 


27.7 


27.6 


0.2 



Table 1: The parameters for the MSSM sample points, MSSMl to MSSM6 in Fig. |3l and the 
breakdown of their contributions to the muon g — 2 vn units of 10~^°. The parameters with 
the mass dimension are given in GeV units, and = is assumed for all the points. The 
numbers in the columns (a) to (e) are the contributions from the corresponding diagrams 
in Fig. [H and the column '(a)-(e)' gives their sum. The numbers in the column 'total' are 
calculated without using the irtz / "^susy expansion, which slightly differ from the sum of the 
five contributions. The last column gives the pull factor, Eq. f|TT]) . 

2.2 The muon g — 2 m. selected SUSY breaking scenarios 

In the previous subsection we have examined SUSY contributions to the muon g — 2 without 
assuming specific SUSY breaking scenarios. In this subsection we examine several SUSY 
breaking scenarios that are consistent with the other constraints like the 6 — ?■ 57 decay rate, 
and discuss in detail their contributions to the muon g — 2. We will later examine their 
predictions for the electroweak observables. 

We take seven scenarios that predict the muon g — 2 values within or very close to the 
1-(T allowed region; a few sample points each from three SUSY breaking scenarios, namely, 
the minimal supergravity (SG) [73], the gauge mediation (CM) [71], and the mixed moduli- 
anomaly (MM) mediation [HHl [73] models. We call those sample points SGI, SG2, GMl, 
GM2, MMl, MM2 and MM3, respectively. 

The SGI point is the mSUGRA sample point advocated as SPSla' in Ref. [76], whose 
main advantage is that it is compatible with all high-energy mass bounds and with the 
constraints from the muon g — 2, Br(6 — 57) and the dark matter relic density. The SG2 
point is a modified version of the SPS4 point, which is a mSUGRA point with tan/3 = 50 
proposed in Ref. [77]. At SPS4 the unified gaugino mass mi/2 is 300 GeV, while at SG2 we 
take mi/2 = 650 GeV so that it is closer to the region favored from the muon g — 2 and 
Br (6 — )■ S7). By this change in mi/2, the pull factors for the muon g — 2 and Br (6 — S7) are 
improved from 3.1 and —5.9 to —0.8 and —1.4, respectively. 

As representatives of the gauge mediation, we take the GMl and GM2 points: in GMl 
tan/3 is large (tan/3 = 42), while in GM2 it is moderate (tan /3 = 15). At these points the 
lightest SUSY particle (LSP) is the gravitino, whose interactions are too weak to be relevant 
for the electroweak observables in the present paper. In the GMl point, which is one of the 
points studied in Ref [78], the next-to-lightest SUSY particle (NLSP) is bino, while in GM2, 
which is suggested as SPS7 in Ref. [77], the NLSP is the stau. Both points fit well with the 
muon g — 2 and Br (6 — )■ 57). 

The MMl and MM2 points are sample points from the mixed moduli-anomaly (MM) 
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tan /3 



Ml 



116 -(88)^ 

465 -(261)2 

214 -(194)2 

120 -(90)2 

255 -(92)2 

108 (103)2 

237 -(102)2 



M2 


Ms 


rriA 


193 


572 


425 


510 


1424 


566 


339 


900 


513 


327 


896 


378 


258 


641 


513 


-248 


-847 


616 


173 


877 


631 



SG 1 
SG 2 
GMl 
GM2 
MMl 
MM2 
MM3 



10 
50 
42 
15 
10 
10 
10 



396 181 

762 585 

504 441 

300 257 

430 188 

-572 253 

534 200 



-445 
-145 
25 
-39 
-465 
245 
509 



103 
277 
181 
169 
170 
-99 
224 



m. 



m. 



m- 



m 



tR {^ihR 



SG 1 

SG 2 
GMl 
GM2 
MMl 
MM2 
MM3 



526 507 505 471 388 

1345 1297 1292 1173 1062 

1329 1269 1263 1264 1165 

861 831 829 836 780 

610 589 546 556 465 

785 796 823 689 585 

758 731 807 705 616 



-(322)2 
-(435)2 

(215)2 
-(241)2 
-(336)2 

(397)2 
-(353)2 



-(153)2 
-(435)2 
-(316)2 
-(153)2 
-(159)2 
(184)2 
-(183)2 



Table 2: The values of the relevant SUSY parameters for the selected scenarios. The param- 
eters with the mass dimension are given in GeV units. (?7ipLR (/ = t, b) are the left-right 

mixing element in the mass-squared matrices of the sfermion /. As for the notation of the 
other SUSY parameters, we use that of Ref . [1] . 





(a) 


(b) 


(c) 


(d) 


(e) 


(a)-(e) 


total 


pull 


SG 1 


25.7 


21.5 


1.5 


-5.2 


-5.4 


38.1 


37.6 


1.4 


SG 2 


20.0 


4.8 


1.0 


-3.4 


-2.8 


19.5 


19.4 


-0.8 


GMl 


34.6 


11.7 


1.4 


-5.3 


-9.2 


33.2 


33.0 


0.9 


GM2 


27.1 


10.6 


1.6 


-5.0 


-9.0 


25.3 


24.8 


-0.1 


MMl 


19.4 


7.2 


1.4 


-4.5 


-1.9 


21.7 


21.7 


-0.5 


MM2 


13.2 


18.8 


0.7 


-2.7 


-4.2 


25.8 


24.7 


-0.1 


MM3 


19.6 


7.9 


1.1 


-3.8 


-1.8 


23.0 


23.1 


-0.3 



Table 3: SUSY contributions to the muon g — 2 for our sample points in units of 10~^°. The 
numbers in the column 'total' are calculated without using the mz/msvsY expansion, which 
slightly differ from the sum of the numbers in the five columns (a) to (e). The last column 
gives the pull factor, Eq. ( ITT]) . 
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mediated SUSY breaking scenario. In MMl, the parameter a, which parametrizes the ratio 
between the moduh and the anomaly mediations, is positive, while it is negative for MM2. In 
MMl and MM2, the parameters (/i, I2, h), which parametrize the contributions from moduli 
to the gaugino masses, are taken to he li = I2 = I3 = I so that it allows the "mirage 
unification" [69], namely the gaugino masses unify at a high scale which can be different 
from the GUT scale ~ 10^^ GeV. In the case of a positive (negative) a, the gaugino masses 
unify below (above) the GUT scale. We take another sample point, which we call MM3, 
from a variant of the MM scenario. At this point, we take {h,l2,h) = (1? 1/2,1) so that 
wino is lighter than binc0. The wino LSP is an interesting possibility since the excess of the 
positron flux observed at PAMELA [82l [83] can be explained by the wino dark matter [S^ . 

For the above seven scenarios we list in Table [2] the values of the relevant SUSY param- 
eters. We only use the parameters in the slepton and ino sectors for the study of the muon 
g — 2, but later we need the squark and the Higgs sectors for the studies of the EW precision 
observables and Br (6 — )■ 57). 

The breakdown of the contributions to the muon g — 2 aX each point with respect to the 
diagrams is given in Table [31 The discussions in the previous section can be verified from 
the numbers in this table. For all our sample points, the diagram (a) gives an important 
contribution. For the points where smuons are relatively light compared to the gauginos or 
the Higgsinos, such as SGI and MM2, the diagram (b) also gives a comparable or larger 
contribution than that of the diagram (a). For all the points, the diagrams (c)-(e) give only 
subdominant contributions. 

The similarity of the SGI,. . . ,MM3 points to MSSMl-6 can be discussed as follows. Since 
SGI is similar to MSSM3 in the sense that it has a bit larger /i than the slepton and the ino 
masses, both diagrams (a) and (b) give important contributions. SG2 is similar to MSSM5 
in /i and the slepton masses but with a heavier inos, and hence the overall size of the SUSY 
contribution is smaller. GMl can be considered to be an interpolation of MS SMS and 6, 
but with a smaller tan/3, and hence the diagram (a) is dominant with a slightly smaller 
contribution from (b). GM2 is a relative of MSSM2, and the breakdown is similar. GMl 
and GM2 have a light right-handed slepton and a moderate-mass (~200 GeV) bino, which 
make the contribution from (e) more important than in other SUSY sample points. MMl, 
MM2 and MM3 are similar to MSSM3, even though they have smaller /i. At MM2, Mi and 
m^^ are smaller than MMl and MM3, which makes the diagram (b) more important than at 
these two points. At MMl and MM3, /i is a bit smaller, and hence the diagram (b) becomes 
a bit less important than at MSSM3 and MM2. 

In summary, similarly to the discussions in the previous subsection, the diagrams (a) 
and/or (b) give important SUSY contributions to the muon g — 2 also in the selected SUSY 
breaking model points. 

* In the original KKLT model [79], the allowed values of (a = 1, • ■ • , 3) are or 1. However, when there 
is a contribution from the dilaton to the gauge kinetic functions, it is possible to have different predictions 
for the gaugino masses from the la = or la = 1 cases [8OI [81] . Here we take into account such a possibility 
by allowing la to take a fractional value as an "effective value" , instead of explicitly introducing the dilaton 
in the gauge kinetic functions. 
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3 The electroweak observables 



In this section, we briefly review the electroweak observables in the framework of Refs. [T||85]. 
and update the parametrizations of the SM predictions. 

The electroweak observables of the Z-pole experiments are expressed in terms of the 
effective Z boson couphngs g^^ [86] to /«, where / denotes the quark/lepton species and a 
stands for their chirahty. The summary of the observables in terms of the effective couplings 
gl^ can be found, for example, in Refs. [HIHT]. A convenient parametrization of the effective 
couplings in generic SU{2)l x f/(l)y electroweak theories is given by [1]: 



where the mean values denote the SM predictions for rrit = 172 GeV, rnHg-^^ = 100 GeV, 
Aa^^^{mz) = 0.0277 and as{mz)5q = 0.118, and the coefficients of Ag^ and As^ control the 
dependences on the oblique (gauge boson propagator) corrections. Here, Ag"^ and As^ are the 
universal gauge-boson-propagator corrections [87] to the effective Z-boson couplings and the 
Z-'-f mixing at the mz scale, respectively, and Ag^. denote the shifts due to vertex corrections. 
In the SM, only {Ag'Dsu and {Ag'}^)sM have non-trivial and niH^^ dependence, and the 
others do not receive rrit- or m//g^- dependent one-loop contribution. On the other hand, all 
the Ag^ terms are non-vanishing at the one-loop level in the MSSM. 

The universal part of the corrections, Ag^ and As^, are defined as the shift in the effective 
couplings and s^(m|) [87] from their SM reference values at {rrit, rriHsM^ ^'^had('^l)) 



The shifts in the two effective couplings can conveniently be expressed in terms of the 
parameters ASz, ATz and x^, 




0.50199 + 0.45250A^| + 0.00469As2 + Agl, 
-0.26920 - 0.24338A^| + 1.00413As2 + Agl, 

0.23207 + 0.20912A^| + 1.00784As2 + Ag%, 

0.34675 + 0.31309A^| - 0.66793As2 + Agl, 
-0.15470 - 0.13942A^| - 0.67184As2 + Ac/]^, 
-0.42434 - 0.38279A^| + 0.33166As2 + Agi, 

0.07734 + 0.06971A^| + 0.33590As2 + Ac/^, 
-0.42116 - 0.38279A^| + 0.33166As2 + Agl, 

0.07742 + 0.06971A^| + 0.33590As2 + A^^, 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 
(12f) 

(12g) 
(12h) 
(12i) 



(172 GeV, 100 GeV, 0.0277): 



^|(m|) = 0.55602 + A^|, 
s\ml) = 0.23048 + As\ 



(13a) 
(13b) 



^fz 
Af 



O.OOAUATz, 

0.00360A5z - O.OO24IAT2 + O.OOOllx, 



(14a) 
(14b) 



Here the parameter 



0.0003 



(15) 
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measures the a{m\) dependence of the effective mixing parameter s 



The parameters 



ATz and ASz denote the shift of Sz and Tz from their values at the SM reference point, 
and are related to the S- and T-parameters as [88] : 

ASz = + ARz, (16a) 
A5, 



a 



/\Tz = AT ^ + 1.49Ai?z. (16b) 

I, and 



The factor 5g is the vertex and box corrections to the muon decay constant, G_f [S 
A(5g is the shift from its SM value, 5g = 0.0055 + A^g [HZ]- The i^^-parameter accounts for 
the difference between S and Sz-, and represents the running effect of the Z boson propagator 



corrections between q 



rrir. 



and = [T]. We define it as 
1 1 



-167r 



(17) 



and /S.Rz denotes the shift from the value of Rz at the SM reference point, 1.1879 [T]: 

i?^ = 1.1879 + Ai?^. (18) 

In this study, we use the IV-boson properties, my/ and Vy/, for the fit. Instead of Af/, as 
the third oblique parameter we take Amj^ = mvi/ — 80.365(GeV) which is given as a function 
of A5, AT, Af/, x„ and A^^, as P 

,A5g 



AmH/(GeV) = -0.288A5 + 0.418AT + 0.337A[/ 



0.0055x„ - 0.126- 



(19) 



We also parametrize the VT-boson decay width, Tw To do so, it is useful to introduce the 
parameter R\y which parametrizes the running of the W boson coupling gwi^f) between the 
zero momentum transfer and = m|, since the decay width is roughly given by 

Tw = 3.3904 X IQ-^mlyGF (l + 8.478 x lO^^Rw + 0.000652;,) 

where, in analogy to Eq. f[T7|) . we define Rw by 

1 1 

and define ARw as the shift from its value at the SM reference point: 

Rw = 2.1940 + ARw. (22) 

The SM contributions to the oblique parameters, Sz, Tz, mw and Rz are given in Refs. [H 
[89] as functions of and rnusM- We update the parametrization as 



R 



w 



— 167r 



(20) 



(21) 



(AS'z)sM 

{ATz)sM 
(Amvy)sM 

(A-Rz)sM 



0.2217X/, - 0.1188x^ + 0.0320x^ - 0.0014x4 + 0.0005a;,, 
-0.0995X;, - 0.2858x^ + 0.1175a;^ + 0.0367x4 + 0.00026^^^ 
-0.0017x,,X4 - 0.0033a;, - O.OOOlxtX,, 
-0.137X/, - 0.019X? + 0.018xj - 0.005a;„ - 0.002a;,, 



(Ai? 



VKjSM 



-0.124 < In 



-0.16 < In 



26 



1 + 



mHsM(GeV) 
23 

m^^g^(GeV) 



-In 








] 









-In 








) 




Viooy 





(23a) 

(23b) 
(23c) 

(23d) 
(23e) 
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T 


T r T 

^EW/QCD [GeV] 


u 


3.1166 


+ 0.0030X, 


3.1377 + 


0.00014xt 


+ 0.0041X, 


0.0000146 


-0.000113 


d, s 


3.1167 


+ 0.0030X, 


3.0956- 


0.00015xt 


+ 0.0019X, 


0.0000032 


-0.000160 


c 


3.1167 


+ 0.0030X, 


3.1369 + 


0.00014xt 


+ 0.0043X, 


0.0000146 


-0.000113 


b 


3.1185 


+ 0.0030a;, 


3.0758 - 


0.00015xt 


+ 0.0028a;s 


0.0000026 


-0.000040 


V 




1 




1 










e,/i 




1 




1 




0.0000368 





r 




1 




0.9977 




0.0000368 






Table 4: The numerical values of the factors C/y, 
expression for the partial widths of the Z boson. 



CfA, ^(jaK ^EW/QCD '^hich appear in the 



The parameters Xt, Xh and Xs are defined as 

mt - 172GeV ln(mj^3^/100GeV) 



3GeV 



In 10 



as{mz)^q - 0.118 
0.003 



(24) 



so that their numerical values are expected to be less than unity. As for the vertex corrections, 
the shift (A^f^^ r)sm in the SM is given by 



(A^^)sM = -0.000058X/, + 0.000128xt, 
(A^^)gM = -0.000042X;, - 0.000025xt 



(25a) 
(25b) 



where the x\ term in the equation for (A5f^)sM is purely from the result of the numerical fit. 

Using the effective coupling (7^, the electroweak observables can be written in the following 
way. First, the Z-boson partial decay width into // is. 



6V2TV 



C 



fv 



fA 




A-^ 

^EW/QCD 



(26) 



The value of each correction factor is summarized in Table |H Cjy and CfA describe the 
corrections to the color factor in the vector and axial- vector currents, respectively, which have 
a dependence on a, and m^. The term S(^^ represents the corrections from the imaginary 
part of loop-induced mixing of the photon and the Z boson. The term A^^^^q^^j-, is the 
non-factorizable mixed electroweak and QCD corrections [9D] , whose values in Table H] have 
been copied from the second paper of Ref. [91]. Qf is the electric charge of the fermion / in 
the normalization that Qf = —1 for the electron. 

As a check of our parametrization, in Fig. H] we give a comparison of Xmm constructed 
from the fit for the first 15 observables in Table [5] together with m^, Aa|j^j(m|), ds(m2)5g 
as a function of rriHsM by using ZFITTER |9T] and that fitted by using our parametrization. 
Since our parametrization is designed so that it gives a good description only in the region 
lOOGeV < rriHs^^ < 1000 GeV, we find that the agreement becomes worse for rriHs^^ ^ 100 
GeV. 

In Table [5l we show the electroweak observables used in the present analysis. The 
experimental values of the Z pole observables, including the correlations among errors that 
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data 


SM 


without Aj 


LEP 1 




best fit 


pull 


best fit 


pull 


line-shape & FB asym.: 












T, (GeV) 
cTO(nb) 


2.4952(23) 
41.540(37) 


2.4958 
41.478 


-0.25 
1.67 


2.4963 
41.478 


-0.46 
1.69 


Ri 


20.767(25) 


20.743 


0.98 


20.746 


0.82 


^FB 


0.01714(95) 


0.01647 


0.71 


0.01690 


0.26 


r polarization: 

A / T "i \ 

Ar{Pr) 


0.1465(32) 


0.1482 


-0.52 


0.1501 


-1.11 


b and c quark results: 














0.21629(66) 


0.21583 


0.70 


0.21581 


0.72 


7-) 

Uc 


0.1721(30) 


0.1722 


-0.05 


0.1723 


-0.05 


^FB 


0.0992(16) 


0.1039 


-2.92 








^FB 


0.0707(35) 


0.0743 


-1.02 








jet charge asymmetry: 












SLC 


0.2324(12) 


0.2314 


0.85 








^LR(^e) 


0.1513(21) 


0.1482 


1.51 


0.1501 


0.58 




0.923(20) 


0.935 


-0.58 










0.670(27) 


0.668 


0.06 








Tevatron + LEP 2 












(GeV) 


80.399(23) 


80.376 


0.99 


80.400 


-0.04 


(GeV) 


2.085(42) 


2.092 


-0.16 


2.093 


-0.19 


Numerical inputs 












(GeV) 
Gf{10 GeV ) 


91.1875(21) 
1.16637(1) 










Parameters 












as(?^z)5q 


0.02759(15) 


0.02761 


-0.14 


0.02759 


0.01 


0.1184(7) 


0.1184 


0.00 


0.1184 


0.06 


(GeV) 


172.0(1.6) 


172.3 


-0.17 


171.9 


0.05 


rriH.., (GeV) 




84.4 




48.9 




Amin 






17.49 




5.35 


d.O.f. 






18-4 




13-4 



Table 5: The electroweak precision data, the SM best fit values and the pull factors. The SM 
predictions and the pull factors have been calculated by using ZFITTER [91] by varying m^, 
mHsMi (^s{fnz)bq and Aa"^l^{m\) as input parameters. The values of the observables are taken 
from Ref. [3], except that niw, Tw, fnt and as{'mz)5q are from Ref. [13], and Aa[2j("^|) is 
from Ref. [8]. The values of and Gp are fixed throughout the calculation. The correlation 
matrix elements of the Z line-shape parameters and those for the heavy-quark parameters 
are found in Ref. [3]. We also show the SM fit and the associated pull factors in the case 
where we do not use the jet asymmetry data, namely, A^^, A^^, sin^ (^^ei^, A^ and Ac- 
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30 50 100 150 200 300 



Figure 4: The comparison of Xmm ^^e electroweak observables as a function of the SM 
Higgs boson mass tuhsm fitted by using our parametrization (sohd hne) and by using the 
output of ZFITTER (dashed hne). Our parametrization is vahd for ttlhsj^ > 100 GeV. 

are not reproduced in Table [5l are taken from Ref. [3]. The values of mw,^w,^t and 
o^s{fnz)5q are taken from Ref. [I3], and the value of Aa^^^{'m'^z) ^ 

^«had("^l) = 0.02759 ± 0.00015, (27) 

is from Ref. [8], in which the prediction Eq. ([2]) for the muon g — 2 is found. In Table [5] we also 
show the SM best fit values calculated by using ZFITTER by varying m^, mnsM^ ^sij^z) and 
Aa^^l^{m\) as the input parameters, mz and Gp are fixed at their central values throughout 
our analysis. 

Here we consider two cases, the case using all data and the case without using the jet 
asymmetry data, namely, A^p^, ^fb' sin^^p*^*, and Ac, because there is stih theoretical 
uncertainty in the calculation of QCD corrections [57]. In the last two rows of Tabled we 
show the values of Xmin and degrees-of-freedom (d.o.f.), which is the number of used data 
minus the number of input parameters. From the fit and the value of Xmin^ can see that 
the SM with the light Higgs boson gives a good description of the data. If we remove the jet 
asymmetry data, a lighter Higgs boson is favored. Once the best fit parameters are fixed, 
the corresponding values for the observables can be calculated immediately, and these SM 
best fit values and the associated pull factors are also shown in Table [5l 



17 



4 The precision data and the MSSM 



In this section, assuming that there is new physics which gives rise to finite corrections to 
/S.Sz and AT^, we estimate the region of AS*^ and AT^ favored by the Z-pole observables. 
Then, under the assumption that the new physics is the MSSM, we use the constraints from 
ASz and AT^ to find the favored range of the MSSM parameters. Later in this section we 
use mw as another observable to constrain the favored SUSY parameters. We conclude this 
section with the discussion of the case where we do not use the jet asymmetry data. 

4.1 Oblique Corrections 

In this subsection, we first identify the favored parameter range of AS^ and AT^. The 
assumptions to compute the theoretical predictions are the following. The input free pa- 
rameters from the new physics are taken to be AS^, AT^, and ^g\. All the other vertex 
corrections ^gl, are neglected for simplicity. As for the SM parameters, we fix Xt and 
at = as a "reference point". Consequences from different choices of Xt and x^ 

can be easily drawn, as discussed later. We take the reference SM Higgs boson mass to be 
rriHsM — 120 GeV in this section. As for the other input SM parameter Xg, instead of fixing 
it at Xs = 0, we include it in the function, and only after finding the minimum of the 
function, we integrate out Xg. 

Using the mean values, the errors and the correlation matrix of the observables in Ref. [3], 
we obtain 



where pcorr is the correlation between ASz and ATz- In Fig. |5l we show the contours for 
the 39% and the 90% confidence levels (CL) as shown in Eqs. f l28|l . and also plot the SM 
prediction for rriHg-^ = 120GeV, = 172GeV and Aa^^^{m'^^) = 0.0277 as the big open 
circle at the origin. Although the ASz and ATz values which give the minimum value are 
slightly different from the prediction at the SM reference point, these shifts are within the 
l-cT error. We also illustrate how the SM reference point moves according to the change of 

from 172GeV to 175GeV and to = 169GeV, as the "ruler" toward the right end of the 
figure. As we can see, the SM prediction for ATz becomes larger for larger rrit, while ASz 
does not change very much because of the stronger dependence of ATz on m^, see Eq. fl23bp . 
The dependence of the plot on Aa^^^{m'^2;) is shown as another "ruler" at the bottom of the 
figure. For example, if Aajj^j(m|) — 0.0277 > 0, then the origin moves to the right, and the 
agreement of the SM reference point to the data becomes better. 

In Fig. m we also show separately the slepton, squark and ino contributions to ASz and 
ATz- In the figure we take tan/3 = 10, but these SUSY contributions do not change very 
much for 10 ^ tan /3 ^ 50. The qualitative behaviors of those contributions on the ASz- ATz 
plane have been studied in Ref. [1]. 

In the figure, the contributions to ASz and ATz from the sum of three generations 
of squarks for the cases mQ=300 GeV and 500 GeV without left-right mixing among the 




0.020 - 1.'2'2Ag\ - 0.031a;„ ± 0.106 
0.053 + 0.50A^^± 0.137 

-0.00033 ± 0.00082, 

15.5, 




0.91 



(28) 



(29) 
(30) 



18 



0.1 



0.08 - 



0.06 - 



-0.02 - 



-0.04 - 



-0.06 - 



-0.08 




1 \ 
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ASz+222Agl 

Figure 5: The squark, slepton and ino contributions to {ASz, ATz) for tan/3 = 10. The 
SUSY breaking scalar masses for the left-handed and right-handed squarks are assumed to be 
same, denoted by uIq. The ti-tn and bi-bn mixings are controlled by Acs = v4*jj = A^^^. The 
left- and right-handed sleptons are also assumed to have a common SUSY breaking scalar 
mass mj;. The reference SM point, (m^, m//g^, Aa[^2i("^l)) = (172GeV, 120GeV, 0.0277), is 
marked by the open circle at the origin of the plot. If a different value of is chosen, then 
the origin would move according to the scale shown at the right-hand side. Similarly, if a 
different value of Aa||^j(m|) is chosen, the origin would move according to the scale shown 
at the bottom. 
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Figure 6: The predictions from our SUSY sample points, superposed on Fig. |5l 
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squarks are given as the filled squares. In the figure, we assume ttiq = rriu^ = mj^ for 
simplicity. The effects of the left-right mixing on these predictions are shown by the dashed 
lines starting from these squares. On each dashed line, the squark soft mass uiq is fixed at 
the same value, and along the dashed line, the parameter Aes which controls the left-right 
mixing is varied from to 300 GeV. (The definition of Acs is the same as in Ref. [T]H ) 
The predictions for A^s = 200 GeV and 300 GeV are shown by the different squares on the 
dashed lines. 

Similarly, the contributions to ASz and ATz from the sum of three generations of sleptons 
in the cases without left-right mixing are given by the filled diamonds labeled as mi = 100, 
120, . . ., 300 GeV. In the figure, for simplicity, we assume mi = m^, similarly to the squark 
case. Attached to these diamonds are the dashed lines which show the effects of the left-right 
mixing on these predictions. On each dashed line, the slepton soft mass mi is fixed at the 
same value, and the size of the left-right mixing is varied by using the lighter stau mass as 
the measure of the left-right mixing. In the figure, the cases where the lighter stau masses 
are 150, 100 and 81.9 GeV are shown by the different diamonds. 

In Fig. [5l the ino contributions are also shown. The filled upward triangles are the 
predictions for the cases where the lighter chargino masses are 110 GeV and 140 GeV, 
respectively, in the Higgsino-like chargino cases with the ratio M2//X fixed at 100. On each 
dashed line the lighter chargino mass is fixed at the same value, while the ratio M2//X is 
varied from 100 to 0.01 along the dashed line. 

The contributions to ASz from the squarks and sleptons can be understood as follows [1]. 
The AS^-parameter is defined as the sum of AS and ARz- When the left-right mixings of 
the sfermions are negligible, to one-loop order, AS* receives contributions from left-handed 
sfermions, and is proportional to the hypercharge Yy of the sfermion / in the loop. The 
sign of the hypercharge is opposite between the left-handed squarks (Yg^ = +1/6) and the 
left-handed sleptons (Y^^ = —1/2), and this determines the sign of ASz in the limit of no 
left-right mixing. On the other hand, the sign of Ai?^-parameter is always negative for both 
squarks and sleptons contributions [T], and it adds up with AS" constructively for sleptons 
while destructively for squarks. This is why ASz is negative for the sleptons and almost 
zero for the squarks. 

The ATz-parameter is also defined as a linear combination of AT and ARz with small 
corrections from 5g — 0.0055. As mentioned above, ARz is negative, but its magnitude for 
the sfermions is tiny compared to the contribution to AT [1]. 

To discuss the contributions to the AT-parameter from the sfermion sector, it is useful 
to separate three cases depending on the size of the left-right mixing of the sfermion: cases 
without the left-right mixing, with small left-right mixing and with large left-right mixing. 
First, in the case without left-right mixing, the contributions from the third generation 

^ For completeness, the definition of AcS is as follows: the left-right mixing elements of the stop and the 
sbottom mass-squared matrices are given by mt{At — /i/tan/3) and mb{Ab — /itan/3), respectively. (We are 
neglecting possible CP-violating phases for ^,At and Ab for simphcity.) We define A*g and A^^f so that 
these left-right mixing elements are equal to mtA^^g and mbA^g, respectively: A^^^ = At — fi/ tan/3, Ajj^j = 
Af, — /itan/3. In this paper, for simplicity wc only consider the case A^^g = A'^f^, and denote this common 
value by A^^g. 
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squarks can be written as 
(AT) 



GpCq i^j-^-J 




(31) 



where Cq is the color factor {Cq = 3 for the squarks) and we take the limit where the 

squarks are heavy compared to m^. The slepton contribution can be obtained by the obvious 

replacements Cq ^ 1, ii ^ u and bi cl- Second, when left-right mixing is small enough, 

the T-parameter decreases as A^s increases, as studied in Ref. [Ij. Third, when in the limit 

that left-right mixing is large, the T-parameter increases as A^g increases [72] . The behavior 

of the stau contribution to ATz interpolates the above two limits. 

The ino contributions are small in general once we impose the experimental constraint 

from the direct searches on the lightest chargino mass, unless the ino masses are close to the 

experimental bounds. When inos are light, the contributions to ARz can be sizable, which 

make negative contributions to Sz and Tz [I]. In Fig. |5l we show the cases m~~ > llOGeV 

for M2/ ^ = 0.01, 1 and 100. For the cases M2/ ^ = 0.01 and 100, the predicted trajectories 

for {ASz, ATz) overlap on the line ATz = 1.49AS'^. This can be understood in the following 

way. Since the wino mass parameter M2 and the Higgsino mass parameter fi do not break the 

SU{2)l X U{1)y symmetry, the contribution from the ino sector to the S- and T-parameters 

can only come from the off-diagonal elements of the ino mass matrices. When M2/ /i = 0.01 

or 100, the mixing of the ino mass matrices are suppressed by m|/(M| — /i^), which is small 

once we impose m~- > llOGeV and the strong hierarchy between M2 and u. Hence for 

these hierarchical cases, the contributions to ASz and AT^ only come from ARz, namely, 

ASz = ARz and AT^ = lA9ARz, which makes the trajectories on Fig. [S] overlap. 

In Fig. [6l we plot our SUSY sample points on the same frame as that of Fig. |5l 

In the cases of the six MSSM sample points, we can ignore the squark contributions 

because we set all the squark masses to be 2 TeV. So, almost all the MSSM scenarios are 

put near the slepton lines or ino lines. Only the sample point 1 is apart from both lines, 

because it has a sizable contribution from A6g to AT^ by 0.036. Although it has the 

slepton contribution at = 300GeV with small mass splitting and the ino contribution at 

M2/ ^ = 0.75 and m~- = llSGeV, these contribution are almost canceled out. The sample 

point 2 has a quite large contribution from slepton, because it has light sleptons. And it shifts 

above the solid line which shows the slepton contribution without the left-right mixing case 

due to the A6g contribution to AT^ by 0.015. The sample points 3 and 5 are characteristic 

in the ino contribution, because of a light chargino mass, m~- ~ 150GeV, and a small M2/ /i 

Xi 

ratio, M2/ jj. = 0.19. In particular, the sample point 5 has essentially only ino contribution 
because of heavy sleptons, ~ SOOGeV. On the other hand, the sample point 3 has a non- 
negligible slepton contribution (m^ ~ 200GeV) compared to the sample point 5, and a small 
contribution from A5g of —0.010 to AT^. The sample point 4 is determined almost only 
by slepton contribution from = 150GeV and rrir^ ~ lOOGeV. The sample point 6 is an 
interesting point, because it is the case of large SUSY breaking mass, rn^ = 300GeV with 
a light stau, m^^ = 144GeV. In this parameter region AT increase as A^g or squared mass 
difference increase. 

As for the predictions from the selected scenario points, for all the selected model points 
except MM3, the ino contributions are negligible because there m~- > 200GeV. The SG2 
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from lept. asym. data only 



from jet asym. data only 



average 



-0.0025 -0.002 -0.0015 -0.001 -0.0005 0.0005 0.001 0.0015 

Figure 7: The favored ranges of As^ from the leptonic and the jet asymmetry data. Also 
shown is the favored value for the case where we combine all the leptonic and the jet asym- 
metry data. 

scenario gives the largest contribution to AT^. This contribution mainly comes from slepton 
with a large with a large squared mass splitting, and the contributions from the other 
sectors are negligible because the squark masses are more than 1 TeV. The GM2 and MM2 
scenarios also have heavy squarks, and in the large SUSY breaking mass with light stau 
region, (m^,mfj = (250GeV, 120GeV) and (245GeV, 104GeV), respectively. The SGI and 
MMl scenarios have similar parameters, and are located at almost the same place in the 
ASz-^Tz plane. At the GMl point, we can neglect the slepton contribution because of 
large mass, ~ 440GeV, and it has the contribution from squarks at rriq ~ SOOGeV and the 
A6g contribution, 0.014. 

In Figs. [5] and O we do not show the contributions from Higgs bosons, which are known 
to be small for ruh > 115 GeV [1]. In fact, for our selected scenarios SGI,. . . ,MM3, the 
contributions from the Higgs sector is < ASz ^ 0.004 and —0.003 ^ ATz < 0. 

4.2 Data without Jet Asymmetry and Oblique Corrections 

It may be worth repeating the above analysis after excluding the jet asymmetry data from 
the input Z-pole precision observables. 

The pull factor in Table |5]shows that the data of fo-jet forward-backward asymmetry differ 
from the theoretical expectation by roughly three standard-deviations. This can be seen more 
clearly if we look into the favored region of As^ separately from leptonic asymmetry and jet 
asymmetry. The results are summarized in Fig. [71 From the figure we can see that the value 
of As^ determined from the leptonic asymmetry data (Apg, A^(P^), A^^^) does not agree 
very well with that determined from the jet asymmetry data (ApB' ^fb' sin^6'|,^*. At,, Ac). 
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In fact, the values for As^ for these cases are separately 



= -0.00035 ± 0.00021 (lepton only), (32) 
As^ = +0.00072 ± 0.00027 (jet only). (33) 

The fit for the lepton asymmetry data gives Xmm/'^-'^-^- — 1-64/(3 — 1), or the probability 
44%. On the other hand, the fit for the jet asymmetry gives 0.39/(5 — 1), or the probability 
98%. The values in Eqs. (15^ and (15^ differ by 3.1 a. If we average the two values blindly, 
we obtain 

As^ = 0.000044 ± 0.00016 (all asymmetry data), (34) 

with Xmin ~ 9-9. This implies that the asymmetry data agree well within the leptonic data 
and the jet data separately, but not very well between the two sets. Although the same result 
f p^ is obtained by averaging all the asymmetry data at once, with Xmm/'^-O-f- — 11-9/ (8 — 1), 
we feel that this low value of Xmin/^-O-f- is an artifact caused by using data with large 
statistical errors. Since we take seriously the possible deviation from the SM in the muon 
g — 2, we would like to take the difference between ( l32i) and ( l33l) seriously. 

Recently, the jet angular distribution in e~^e~ annihilation has been re-examined [57] in 
the framework of soft-collinear effective theory [92] and a local current-three-parton (qqg) 
operator which contributes to the reduction of the forward-backward asymmetry has been 
identified, and the associated parton shower (jet function) has been obtained in the NLL 
approximation of massless QCD. Although the quantitative effect estimated in Ref. [57] 
reduces the discrepancy between the quark and lepton measurements only slightly, the ob- 
servation suggests that we may need to develop a parton shower program which is capable 
of simulating the jet angular distribution with the accuracy matching that of the precision 
measurements. Until the data can be re-analyzed by using such advanced tools, it may be 
worthwhile to examine consequences of dropping the constraints from all the jet asymmetry 
measurements. 

When we leave out the data for A^^, A^ps^ sin^^^J^*, Ai, and the favored region 
becomes 



ASz = -0.109 + l.bOAgl - OMlx^ ± 0.113 
ATz= 0.024+ 1.24A^i± 0.137 



0.87, (35) 



and the value of the minimum of is 



^'"''^ V 0.00083 J ^ ^ 

As shown in Fig. |S1 we find that the favored region has been shifted to the negative ASz 
direction. We also overlay the MSSM predictions already discussed in Fig. O We find that 
the favored region can be reached by relatively light (~ 100 — 120GeV) sleptons. It is 
interesting to note that these light sleptons can explain the muon g — 2 anomaly naturally. 
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Figure 8: The same figure as Fig. |6] except that we omitted the jet asymmetry data. 
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80.354 GeV is the SM prediction for the SM reference point, [m^^mH, 
(172GeV, 120GeV, 0.0277). The dependences of the SM prediction on and Aa 
are shown at the bottom. Also shown as the red/solid line and the yellow/shaded band are 
the experimental mean and the l-a uncertainty, respectively. 
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4.3 mw in the MSSM 



In our framework, the W boson mass is a quantity which can be calculated from input param- 
eters. The predicted SM value of mw for our SM reference point (m^, m^^g^, Att[|^^j(m|)) = 
(172GeV, 120GeV, 0.0277) is given in Fig. M as the vertical dashed line at mw = 80.354 
GeV. We see that it is away from the experimental result, whose mean is shown as the solid 
vertical line at mw = 80.399 GeV together with its uncertainty shown as the band, roughly 
at 2-a level. At the bottom of the figure we show the dependences of the SM prediction 
for mw on mt and AQ;[2^(m|). When m^ becomes larger the prediction also becomes larger 
because Amw has a rather strong dependence on AT, see Eq. f|T9|) . and AT has also positive 
dependence on m^ as Eq. (123bp . The dependence on Aa[^^|j(m|) is not as strong as on mt, 
but is not negligible. 

In Fig. [9l we also show the individual contributions to the W boson mass from each 
sector in the MSSM. In the figure we take tan /3 = 10, but these SUSY contributions do 
not change very much for 10 ^ tan 50. The squark and slepton contributions to mw 
are examined for the same parameter space in Fig. [61 They make the fit to the mw data 
better than the SM. The sfermion contributions shift mw into the l-cx favored range, when 
rrig <, 500GeV or when mj^ ^ 140GeV. These improvements mainly come from the AT and 
the AS* terms in Eq. f[T^ for the squark and the slepton, respectively. As for the dependence 
on the Acs term, in the case of squarks, a larger left-right mixing makes the correction to 
mw smaller. This can also be explained by the dependence of AT on A^s, as already seen in 
Fig. [6l Similarly, also for the sleptons, when the left-right mixing is not extremely large, the 
larger A^s predicts smaller Amw- However, when the left-right mixing is extremely large, 
the contribution to AT becomes large, as already discussed, which also makes Amw large, 
as seen for (m^, J = (300, ~ 100) GeV in Fig. [H 

The ino contributions to mw are examined for llOGeV < m~- < 140GeV and for 
M2/ jj. = 0.01, 1 and 100. They are relatively small compared to the squarks and sleptons. 
Among the three cases, only the mixed case {M2/^ = 1) gives a sizable correction to mw- 
This can also be understood from the discussion on AS and AT in Section 4.1. 

In Fig. [9] we also show the predictions from the sample SUSY parameter sets. We see 
that for all the SUSY sample points the predicted values for mw are improved compared 
to the SM reference point. Among them, when there is a light slepton, like at MSSM2, 
MSSM4 and MSSM6, the improvement is large since the slepton contributions are larger 
than inos for similar masses. In particular, at MSSM2, where both the sleptons are light, 
the improvement is most effective. 

Also for the predictions from the selected SUSY breaking scenarios, SGI,. . . ,MM3, the 
points with light sleptons make large contributions to mw, like SGI, GM2, MMl-3. At SG2, 
the left-right mixing of the stau makes the contribution large. At GMl, since there are no 
light sleptons, the contribution is small. 

In Fig. |9] we do not show the dependence on the Higgs sector since it is known to be 
small [T]. For our selected scenarios SGI,. . . ,MM3, the contributions from the Higgs sector 
is — 3MeV^ Amw < 0, which is negligible compared to the experimental uncertainty. 
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4.4 Fw in the MSSM 



The SUSY corrections to the W boson decay width, Tyy, can also be calculated once the 
SUSY parameters are fixed. The SM prediction for {mH^^.mt) = (120, 172) GeV is F^y = 
2.090 GeV, which is consistent with the experimental value, = 2.085 ± 0.042 GeV [T3] . 
Compared to the experimental uncertainty, the SUSY corrections to T^^ are very small 
(~ 0.001 GeV) for our sample SUSY parameters, and we find that F^/ is not as useful as 
other EW precision parameters to constrain SUSY contributions. 

4.5 Summary of Electroweak Observables 

In Tables IH] and [7] we summarize the values of the SUSY contribution to the oblique param- 
eters and the electroweak observables for our sample parameters. From the tables, we see 
that for our sample parameters the SUSY corrections are small in general since for those 
points the SUSY particles are at the range of a few hundred GeV or heavier. 

Next, if we look into the observables in Tables [6] and [TJ we see that the observables like Ac 
and sin^ O]!^^^ do not depend on SUSY parameters very much compared to the experimental 
accuracy. We also see that some jet asymmetry observables like Ap^, Ap^ and Af, do not 
agree between experiment and SM, which SUSY contributions cannot improve very much, 
as is well known. 

We also give Xmin ^^e cases where (i) all the data are used, (ii) only Ap^ is excluded, 
(iii) the jet asymmetry data (data with * in the tables) are excluded. We see that the values 
of Xmin show sizable changes between the cases (i) and (ii), but not very much between (ii) 
and (iii). This may suggest that A^^ is the main source of the deviation of the SM from the 
data. 

In Tables E] and [7] we also give the SUSY contribution to the shift Arriw. Since the shift 
can be written in terms of the oblique parameters by Eq. f|T9|l , we can calculate the shift also 
from the values of AS*, AT, and so on in Tables [6] and [3 Since for our sample points AT 
is larger than AS, AU and — A^c/a, AT gives the main contribution to Amiy. In those 
sample points with larger AT, such as SGI and SG2, the predicted shift Amw is larger, 
which makes the fit of rriw better. Similarly, for those point with smaller AT, such as GMl 
and MM2, the shift Amw is small, which makes the total worse. 

The SUSY contribution to Tw is small compared to the experimental accuracy, as seen 
from Tables [6] and [3 We conclude that it is not very useful to constrain SUSY contributions. 

In this paper, we do not consider the SUSY non-oblique corrections other than ASq- In 
the MSSM we expect that the corrections to the Z-b-b vertex is the largest among vertex 
corrections. We find that the SUSY contributions to Ag\^i^ are at most of the order of 10~^, 
which is far smaller compared to the oblique corrections. 
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data 


SM 


MSSMl 


MSSM2 


MSSM3 


MSSM4 


MSSM5 


MSSM6 


AS 






0.033 


-0.025 


-0.001 


-0.010 


0.009 


0.029 


AT 






0.040 


0.048 


0.016 


0.023 


0.005 


0.038 


ARz 






-0.032 


-0.017 


-0.023 


-0.007 


-0.020 


-0.002 


Amw (GeV) 






0.017 


0.032 


0.009 


0.013 


0.002 


0.009 


Tz (GeV) 


2.4952(23) 


2.4948 


2.4954 


2.4960 


2.4945 


2.4952 


2.4943 


2.4950 


< (nb) 


41.540(37) 


41.481 


41.494 


41.477 


41.488 


41.481 


41.486 


41.482 


Ri 


20.767(25) 


20.737 


20.734 


20.745 


20.731 


20.739 


20.733 


20.735 


n / 


0.01714(95) 


0.01613 


0.01626 


0.01651 


0.01622 


0.01628 


0.01614 


0.01611 




0.21629(66) 


0.21585 


0.21587 


0.21586 


0.21586 


0.21585 


0.21578 


0.21578 


Rc 


0.1721(30) 


0.1722 


0.1722 


0.1722 


0.1722 


0.1722 


0.1722 


0.1722 


^f'b 


0.0992(16) 


0.1028 


0.1032 


0.104 


0.1031 


0.1033 


0.1029 


0.1028 




0.0707(35) 


0.0735 


0.0738 


0.0744 


0.0736 


0.0738 


0.0734 


0.0734 




0.2324(12) 


0.2316 


0.2315 


0.2314 


0.2315 


0.2315 


0.2316 


0.2316 


* ^fe 


0.923(20) 


0.935 


0.9347 


0.9348 


0.9346 


0.9347 


0.9354 


0.9354 


* Ac 


0.670(27) 


0.668 


0.668 


0.669 


0.668 


0.668 


0.668 


0.668 


AriPr) 


0.1465(32) 


0.1467 


0.1473 


0.1483 


0.1471 


0.1473 


0.1467 


0.1466 


^LR(^e) 


0.1513(21) 


0.1467 


0.1473 


0.1483 


0.1471 


0.1473 


0.1467 


0.1466 


mw (GeV) 


80.399(23) 


80.354 


80.370 


80.386 


80.363 


80.367 


80.356 


80.362 


TvK (GeV) 


2.085(42) 


2.090 


2.092 


2.093 


2.091 


2.091 


2.091 


2.091 


xlw (all) 




20.61 


18.16 


18.37 


20.07 


19.21 


21.74 


20.39 


Xew (excl. Apg) 
Xlw (excl. *) 




16.32 


12.76 


10.59 


15.08 


13.71 


17.26 


16.11 




14.83 


11.02 


8.30 


13.43 


11.94 


15.72 


14.61 



Table 6: The breakdown of the radiative corrections at our SUSY sample points MSSMl - MSSM 6. In the column labeled as "SM", the 
values for the SM reference point (m^, m//g^, Aa[|3^j(mD) = (172GeV, 120GeV, 0.0277) are given. The values for the SUSY corrections 
A5, AT, ARz and Am\Y are the deviations from this SM reference point. 





SGI 


SG2 


GMl 


GM2 


MMl 


MM2 


MM3 




0.008 


0.012 


0.009 


0.006 


0.003 


0.005 


0.003 


AT 


0.044 


0.054 
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0.034 


0.019 


0.027 


Amw (GeV) 


-0.019 


-0.002 


-0.005 


-0.009 


-0.011 


-0.011 


-0.019 


0.019 


0.020 


0.003 


0.013 


0.015 


0.008 


0.013 


T, (GeV) 


2.4958 


2.4957 


2.4948 


2.4955 


2.4956 


2.4952 


2.4952 


< (nb) 


41.486 


41.484 


41.485 


41.483 


41.483 


41.485 


41.486 


Ri 


20.738 


20.734 


20.734 


20.739 


20.739 


20.736 


20.736 


, n / 

^FB 


0.01631 


0.01630 


0.01615 


0.01626 


0.01629 


0.01620 


0.01627 




0.21589 


0.21569 


0.21575 


0.21588 


0.21588 


0.21587 


0.21588 


Rc 


0.1722 


0.1723 


0.1722 


0.1722 


0.1722 


0.1722 


0.1722 


* ^FB 


0.1034 


0.1035 


0.1030 


0.1032 


0.1033 


0.1030 


0.1032 




0.0739 


0.0739 


0.0735 


0.0738 


0.0738 


0.0736 


0.0738 


* sin^ 9^^l^ 


0.2315 


0.2315 


0.2316 


0.2315 


0.2315 


0.2315 


0.2315 


* ^6 


0.935 


0.936 


0.936 


0.935 


0.935 


0.935 


0.935 


* 


0.668 


0.668 


0.668 


0.668 


0.668 


0.668 


0.668 


Ar{Pr) 


0.1475 


0.1474 


0.1468 


0.1472 


0.1474 


0.1470 


0.1473 


^LR(^e) 


0.1475 


0.1474 


0.1468 


0.1472 


0.1474 


0.1470 


0.1473 


vfiy^ (GeV) 


80.373 


80.373 


80.357 


80.366 


80.369 


80.362 


80.366 


(GeV) 


2.092 


2.091 


2.090 


2.091 


2.091 


2.091 


2.092 


xlw (all) 


18.08 


19.38 


21.44 


18.97 


18.65 


19.70 


19.05 


Xew (excl. Apg) 
xlw (excl. *) 


12.19 


13.23 


16.69 


13.53 


13.01 


14.85 


13.61 


10.34 


11.30 


15.08 


11.78 


11.22 


13.24 


11.86 



Table 7: The breakdown of the radiative corrections at our sample points in the selected SUSY models. 
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Figure 10: The favored region in the plane of the most experimentally constrained linear 
combination of ASz and ATz versus myy. The inner/outer ellipses stand for the 39%/90% 
CL favored regions. The upper open circle is the SM prediction for (m^, rriHs^, A«[^^^(m|)) = 
(172GeV, 120GeV, 0.0277). Also shown as the red/filled blobs are the predictions of our 
sample SUSY models. The lines with ticks around the lower open circle are the "rulers" to 
show how the SM prediction, and hence all the SUSY model prediction points as well, shifts 
when more accurate data on rrit and Aajj^|j(m|) are available. 

5 Preferred parameters in a few SUSY breaking mod- 
els 

In the previous section, we have studied the constraints from the EW precision data on 
ASz, ATz, and m^y. In Figs. E] and [HI the constraints on ASz and ATz are shown by 
the favored region of an elliptic shape with a strong positive correlation. This means that 
the linear combination of ASz and ATz along the minor axis of the ellipse is constrained 
much stronger than the orthogonal combination along the major axis. We also note that 
this combination along the minor axis direction is strongly affected by the removal of the 
jet asymmetry constraints. In this section, we show the constraints on our g — 2 favored 
SUSY models from the electroweak precision measurements in the plane of this strongly 
constrained combination and rriw, as two-dimensional 'summary plots'. For those models of 
SUSY breaking where the squark masses are related to the slepton and ino particle masses, 
we also examine the constraints in the plane of muon g — 2 and Br (6 — 57). The MSSM 
model points do not appear in those plots since squark masses can be set large to make them 
consistent with the Br (6 57) constraint. 

In the analysis we have performed in Fig. El the direction of the minor axis of the ellipse, 
which corresponds to the most tightly constrained direction, is — 0.799AS'^ + 0.602ATz, 
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Figure 11: The same as Fig. [TO] except that the jet asjTiimetry data are not included in the 
analysis. 

along which 

-0.799A^z + 0.602AT^ = 0.016 + 2.07Ac/^ + 0.025x„ ± 0.036. (37) 

We can conveniently combine this result with the constraint from mw in a single figure 
(Fig. [To]) . In the figure, the 39% and the 90% CL favored regions are shown as the el- 
lipses. Also shown as the upper open circle is the SM prediction for our reference point, 
(mt,mHsM, Aaiad(^l)) = (172GeV, 120GeV, 0.0277). The SM predictions for different m^ 
within 169GeV < rrit < 175GeV can be read off using the "ruler" around the lower open 
circle. For example, if we take = 175GeV instead of = 172GeV, the SM prediction 
moves from the upper open circle toward the upper- right, in the direction of the vector whose 
initial point is the lower open circle with the terminal point being the point shown as "175" 
and also by the length of the same vector. We see that, within the range of the top quark 
mass shown, a larger top quark mass is favored from the data. This preference for a larger 
top quark mass is clearer in rriw than in the most constrained combination of ASz-^Tz- 

In Fig. [TUl we also plot the sample SUSY points. In the direction of the most constrained 
AS'^-AT^ combination the fit does not improve very much by introducing SUSY particles 
since the SM already gives a good description. It is also seen that all our SUSY sample 
points lie within the 1-a favored range of this most constrained direction. These tendencies 
could already be read off in Fig. [61 In the direction, we have improvements in general 
as discussed in Section 4.3. 

We can repeat the same analysis as Fig. [TD] also in the case where the jet asymmetry 
data are not included in the input data. The most constrained direction in the A^^-AT^ 
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plane in this case is 

-0.780ASz + 0.626ATz = 0.100 - 0.39Agl + 0.024x„ ± 0.044. (38) 

Combining this with the constraint from rriw, we show the favored region in Fig. [Ill Com- 
pared to Fig. Uni we see that the eUipses move toward the right, which is the negative ASz 
direction. We also show the SM prediction for the reference point as the upper open circle, 
which will move according to the "ruler" around the lower open circle for a different m^. 
In this case, there is a clearer tendency that the large top quark mass is favored within 
the range of rrit shown. We also show the predictions from our sample SUSY points. As a 
general tendency, our SUSY scenarios slightly improve the fit over the SM reference point. 
In this case, the light slepton scenario, MSSM2, can improve the fit most efficiently among 
our sample SUSY points, since we have chosen the slepton mass small in such a way that it 
can better explain the negative ASz- The degree of the improvement can also be seen in the 
without jet asymmetry data of Tables |6] and [71 At MSSM2, Xmin = 8.30, which is much 
better than that of the SM reference point, Xmin = 14.83, and also than other sample SUSY 
points. 

Having discussed the EW constraints, we now take the constraint from 6 — )■ 57 into 
account. The experimental value quoted in RPP 2010 [13] is Bt{B — )■ Xs'~f) = (3.55 ± 
0.26) X 10-^ while the SM prediction at NNLO is Br(6 ^ 57) = (3.15 ± 0.23) x 10"^ in 
Ref. P], and Br(6 S7) = (2.98 ± 0.26) x 10"^ in Ref. Since the experimental and 
the theoretical values agree within 2-a level, it is preferred that the SUSY contribution is 
not very large so that it does not spoil the rough agreement. 

To suppress the SUSY contributions, we can think of two possibilities: either the relevant 
SUSY particles are heavy enough, or a cancellation among the relevant diagrams happens. 
To see how this can be realized, let us look into the structure of the SUSY contributions. At 
one-loop, the SUSY contribution mainly comes from the chargino-stop loop and the charged- 
Higgs-top loop diagrams. We neglect the possible contributions from the gluino-sbottom 
loop diagram, assuming the minimal flavor violation [95]. Under this assumption, we are 
only interested in those parameter sets in which the chargino-stop and the charged-Higgs- 
top contributions cancel with each other to the extent that the experimental constraint is 
satisfied, or those parameter sets in which the relevant SUSY particles are heavy enough. 

In Fig. [121 we show the predictions for the muon g — 2 and Br(6 — 57) for these sample 
points. Also shown as the inner and the outer ellipses are the 39% and the 90% CL contours, 
respectively. In the figure, to specify the favored region of Br (6 — )■ 57), we use the experi- 
mental result Br(5 — )■ ^^7) = (3.55 ± 0.26) x 10~^. As for the uncertainty in the Standard 
Model prediction, we assign 0.26 x 10"^, which is the larger of the uncertainties in the two 
SM predictions mentions above. We add the uncertainties in the experimental results and 
the SM prediction in quadrature. Concerning the MSSM prediction for Br (6 — > 57), we use 
micromegas version 2.0.7 [96], while the SUSY contribution to the muon g — 2 is calculated 
by using our own code. 

From the figure, concerning Br (6 — > 57), we see that all the points are within the 90% CL 
favored region. We also see that Br(6 — 57) at SGI is a little bit farther from the central 
point than the other points. This happens since, at SGI, the third generation squarks 
and the lighter charginos are slightly lighter than those of the other points, and since the 
cancellation among the SUSY diagrams are milder. 
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Figure 12: The branching ratio of 6 — >■ 57 and the SUSY contribution to the muon g — 2 
at our sample points. The inner and the outer eUipses stand for the 39% and the 90% CL 
favored regions, respectively. 

Concerning the muon g — 2 ioi those sample points, since we have already discussed in 
Section 2, we do not repeat it here. 

As an additional constraint on these SUSY sample points, we now comment on the dark 
matter relic density predicted from these models. At all our SUSY sample parameters, the 
lightest SUSY particle (LSP) is stable because of the -R-parity conservation, and hence the 
LSP is a potential candidate for dark matter. The LSP is the lightest neutralino in our 
sample points based on mSUGRA or the mirage mediation (MM), while those based on the 
gauge-mediated models the LSP is gravitino. In either cases, the relic density of the LSP 
is calculable. For the mSUGRA and MM based points, we have calculated the LSP relic 
density using micromegas. The results are Qh"^ = 0.08, 0.01, 0.11, 0.08 and 0.001 for SGI, 
SG2, MMl, MM2 and MM3, respectively. For all these cases except MM3, the LSP is an 
almost pure bino, with a very small mixture from Higgsinos and wino. For MM3, since the 
LSP is wino, the relic density is smaller because of the larger annihilation cross section of 
the wino LSP. As for the gauge mediated model sample points, GMl and GM2, the LSP 
mass, namely the gravitino mass, is in the eV range, in which case the LSP relic density 
is negligible. These relic densities should be compared to the results of a global fit of the 
cosmological parameters on the non-baryonic matter density Qj^uh"^ |13j . 

^DM/i^ = 0.110 ±0.006. (39) 

We see that for all our sample points, the relic density of the LSP is nearly equal to or less 
than the observed density of dark matter. If the relic density of the LSP were significantly 
larger than the value in Eq. (139|) . those models would be excluded. On the contrary, if 
the relic density of the LSP is less than the value in Eq. (1391) . such a model can still be 
phenomenologically viable since there is still a possibility that an unknown particle like an 
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axion can also contribute to the dark matter density. Hence we conclude that our sample 
parameters are not excluded from the dark matter density calculations. 

We do not include constraints from the low-Q^ precision measurements such as atomic 
parity violation and neutrino-nucleon scatterings at low energies since constraints from these 
measurements are known to be much less stringent than those from the Z-pole experiments. 
Another class of observables we do not include is those from S-physics, namely Bt{B^ 

Bt{B~^ — )■ T'^h'r) and Am^. Even though these observables potentially give non- 
trivial constraints on large tan P models [97], we do not include them since our main interest 
in the present paper is in the signal from the slepton and the ino sectors rather than the 
squark and the Higgs sectors. 

6 Summary 

We have studied impacts of recent muon g — 2 measurements and the LEP final electroweak 
data on the MSSM. We identify several regions of the MSSM parameter space which fill 
the gap between the SM prediction and the observed value of the muon g — 2, and at the 
same time have observable effects for the electroweak precision measurements. In all the 
selected regions of the MSSM parameter space, the MSSM predictions are consistent with 
the LEP/SLC Z boson observables, while improve the SM fit to the W boson mass slightly. 
When we remove the constraints from the jet asymmetry measurements at LEP/SLC, we 
find that MSSM models with very light sleptons 200 GeV) and moderately heavy ino 
particles (~ several 100 GeV) are favored over models with a very light chargino (~ 100 
GeV) and moderately light sleptons (~ a few 100 GeV). 

We also examined a few models of SUSY breaking scenarios, including minimal SUGRA 
models, gauge mediation models, and the mixed moduli and anomaly mediation models. 
All of them have parameter region with relatively light slepton and ino particles which 
contribute to the muon g — 2. Those models with moderately heavy smuons and ino particles 
can still contribute to the muon g — 2 with large tan /3 40), and can at the same time 
improve the fit to mw and the Z boson parameters if there is a significant mixing in the 
stau sector. Sample scenarios in each SUSY breaking models are found which improves 
the SM fit to the muon g — 2, rriw, the Z parameters in the absence of jet asymmetry 
data, and are still compatible with Br (6 — t- 57). We believe that our investigations will 
help us identifying the supersymmetry breaking scenario once signatures of SUSY particle 
productions are discovered at the LHC. 
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